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ABSTRACT 


The fluctuations in the intensity of radio source PKS 2025-15 
during its occultation by Comet Kohoutek (1973f) on January 5, 1974, 
can be interpreted as scintillations due to the turbulent plasma 
in the Comet's tail. It is found that the rms fluctuation of the 

_3 

electron density in the Comet's tail, AN, is about 200 electrons cm , 
the inner scale of the fluctuation is about 8 x 10^ cm and the 
largest scale of fluctuation may reach 6 x 10^ cm. 
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I. INTRODUCTION 


Comets are known to develop plasma tails as they approach the 

Sun. Recently Ananthakrishnan et al. (1975) have reported the 

fluctuations in the intensity of the extragalactic radio source PKS 

2025-15 observed during its occultation by the tail of Comet Kohoutek 

(1973f) on January 5, 1974. The observations were made at 327 MHz. 

They found that the scintillation index M z (defined as( variance 
1 /2 

of intensity) /mean intensity) was ^20% and reached a maximum of 
<30%. The average observed time scale of fluctuations is about 
10 second. However, they concluded that it is difficult to explain 
the fluctuations in terms of scintillations produced by the passage 
of radio waves through the cometary plasma. 

The purpose of this note is to point out that the fluctuations 
of radio intensity observed during occultation can be explained by 
scintillations caused by the cometary tail if we use the scintillation 
theory correctly. In Section II, we use a Gaussian spectrum for the 
plasma irregularities in the cometary tail to analyze the scintillation 
data, while in Section III, a mure realistic spectrum, the Kolmogorov 
spectrum, is used. 
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II. Scintillations During Occultation 


Since the size of plasma tail of a comet is much smaller than 
the distance between the comet and the observer on the earth, we 
will use "thin phase-changing screen approximation" (Salpeter 1967) 
to analyze the problem. A schematic sketch of this approximation is 
presented in Figure 1, in which z is the direction of wave propagation 
and p = (x,y) is the transverse coordinate. The random phase 

a. 

fluctuations of the wave are produced by a layer of irregularity in 
the plasma tail. After passing the screen, as shown in Figure 1, 
the phase of the wave is randomized and is characterized by the 
function <j>(x,y). The intensity fluctuations 6l(z,p) of the 
electromagnetic wave at position z > o are then built up by diffraction 
as the wave propagates to the observer. 

Let the two-point correlation function of phase be 
$2 P ( p ) = <<|)(p 1 ) $ (p t + P )> (1) 

0 O. % % »V> 

with P a (o) = 1. Here < > denotes an average over an ensemble 

<p 

of random phases. PJp) is the normalized phase correlation function 


and <j> 0 is the rms phase fluctuation. Let the normalized intensity 
correlation function of the observed wave be 


<61(0!) 6l (p j + p)> 

Pr ( p ) = -1 


I 




( 2 ) 


where < T > is the mean intensity and 61 is the fluctuation part. The 
scintillation index M z measured at z is defined as 

M z 2 = Pj(o) (3) 
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Now we consider the case in which the phase correlation function 
is Gaussian with a characteristic scale a, i. e. 

P (p) = exp (-p 2 /2a 2 ) (4) 

9 'V, 

Let the correlation scale of the observed intensity fluctuation be 

A A 

Z and let P A (q) and P T (z,q) be, respectively, the Fourier transforms 
of PJp) and P T (z,p) in transverse coordinates. 

9 'V L \ 

The results of the "thin phase-changing screen approximation" 
can be summarized as (Salpeter 1967, Lee and Jokipii 1975b). 

Case (i ), for <t> Q < 1 

P, (2,n) = 4 * 0 2 sin 2 (fj^) (q), l c = a; (5) 

Case (i i ) , for <p Q >_ 1 

(a) P, (z,q) = 4> 2 (z q 2 /k) 2 P (q), Z = a, for z « £ f (6a) 

(b) P T (z,p) = exp [-24 . 2 (I-P.(p))] - exp [-2* 2 ], Z = a/6, for z »Z f . (6b) 

i % 0 <J) ^ UCO T 

Here Z f = ka 2 /<t> Q is the "focal length" and z is the distance from the 

phase screen. For the region where z 'v Z f , the correlation scale Z c changes 

from a to a /<t> o smoothly (Lee 1974). 

For the occultation of radio source PKS 2025-15 by Comet Kohoutek 

(1973f) (Ananthakrishnan et al. 1975), we have: z = 1.3 x 10 cm, M z = 

0.2 - 0.3, and the average time scale of intensity fluctuation (t ) is 

'v- 10 sec. With the transverse velocity of the plasma tail to the line 

of sight, V-j ^ 80 km/sec, we have the spatial correlation scale £ c = 

800 km. Our problem does not belong to Case (i), where 4> o _< 1, since 

[z/(kf 2 )1 = 0.033 would give us M < 0.1. The observed value of M is 
c z z 

between 0.2 and 0.3. It cannot be Case (iib), because z << k Z^ 2 < Z 
We will see that Case (iia) indeed gives us a consistent explanation. 
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D 


For <(> >_ 1 and z << £ f , we have t - a = 800 km. By Eq. (6a), we 

0 " I V 

have 

"l* * V <k^> 2 < 7 > 

With the observed values of z/ka 2 = 0.033 and M z = 0.2 - 0.3 

Eq. (7) gives us 

<f> =2-3 (8) 

o 

Note that with the value of <t> in Eq. (3), L - 2 x 10 14 cm and 

0 T 

z << l* is satisfied. Thus for a = 800 km, and <j> = 2 - 3, we can 
t o 

indeed obtain a scintillation index M of values between 20% and 30%. 

The rms phase fluctuation $ can be related to the rms electron 
density fluctuation, AN, and the thickness of plasma layer, D. We 
have (Lee and Jokipii 1975b) 

<p 2 = 8X2 1 / 2 TT 5 / 2 a r 2 k" 2 D AN 2 (9) 

o e. 

where r is the classical electron radius and k is the wave number, 
e 

Puttinn $ = 2.5, we obtain from Eq. (9) 

D AN 2 = 2.4 x 10 17 cm" 5 (10) 

From Figure (2) of Ananthakrishnan et al. (1975), D is about 6 x 10 10 
cm. Thus we have 

AN = 20 cm' 3 . (11 ) 

The above result must be considered as the lower bound of the rms 
electron density fluctuation since for the Gaussian spectrum, the 
value of AN corresponds only to the density fluctuation of those in- 
homogenei ties which are important in causing scintillations. 

The Gaussian spectrum is rather artificial in interpreting the 
scintillation data and the power-law spectrum is a more realistic form. 
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In the earth's atmosphere (Tatarskii 1971), in the solar wind (Jokipii 
1973), or in the interstellar medium (Lee and Jokipii 1976), the 
spectrum of the turbulent medium has been observed to be close to the 
Kolmogorov power-law spectrum. In the next section, we use the Kolmogorov 
spectrum to analyze the scintillation data. 
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III. Kolmoqorof Spectrum 


Consider the case in which the turbulent cometary tail has a 
Kolmogorov spectrum. It can then be shown that the phase correlation 
function also possesses a Kolgomorov spectrum (Lee and Jokipii 1975a). 
Thus we write 


P,(q) ■ A (Hq 2 L 2 )' a/2 exp(-qV/2) (12) 

with L »l. This spectrum is flat for q < L _1 , is a power-law 
with index - a for L’"* < q < and is cut off for q > i”\ 

Here L is the coherence scale (or the largest scale) and i is termed 
the inner scale. Usually 3 < a < 4. In this section, we set a = 11/3 
corresponding to the Kolmogorov spectrum. In Eq . (12), A is a constant 
and is determined by the condition that P. (0) * 1. For L » i , 
we have 

A * ( | - 1 ) L 2 / n . (13) 

o 

Since Pj(0) = M z = 0.04 - 0.09 from the occultation data, the 
scintillation is still weak and we have 


Pj( z , q) M / sin 2 ( zq 2 /2k ) P (q ) (14) 

(Jokipii 1970, Lee and Jokipii 1975b). It can be shown from Eq.(14) 
that the correlation scale of Pj is 


t 

c 


r i , for t > (z/k)* 2 

- 

i. (z/k)**, for t < (z/k) 1 *. 
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u 

(' 


Noting i f (z/k) In our case, we conclude that i , • i . Thus 
c c 

the inner scale of the turbulence in the cometary tail is 8 x 10^ cm. 

Again the rms phase fluctuation <j> Q is related to aN and the 
thickness D (Lee and Jokipii 1975b) by 

4 > 2 = 76.8 r 2 k" 2 D AN 2 r(l 1/6)/ r(l/3). (16) 

From Eq. (14), we have 

M 2 = Pj (0) = A <> 2 (z/k) 2 8 w r(7/6) f 7/3 L _11/3 (17) 

cy noting that (z/k) ls « i . Assuming that the coherence scale is of 
order of the width of the cometary tail and setting L = D * 6 x lO^cm 
we have from Eq. (1 7) 

AN = 200 cm" 3 (18) 

for M z = 0.25. Note that from Eqs. (13), (16) and (17), aN 2 « i 7/3 L 2/3 

The value of aW in Eq . (18) can be considered as the upper bound for 
the rms electron density fluctuation since L can not be larger than 
the size of the cometary tail. 
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IV. Conclusion 


From the above analysis, we find that for a Gaussian spectrum, 
the rms electron density fluctuation AN of the plasma irregularities 

_3 

in the cometary tail is about 20 cm and the correlation scale a is 

7 -3 

about 8 x 10 cm. However for the Kolmogorov spectrum, aN * 200 cm , 

the inner scale of the turbulence £ = 8 x lO^cm, and the coherence 

scale L = 6 x lO^cm. Since the Kolmogorov spectrum is a more realistic 

_3 

one, we believe that AN = 200 cm is the better value for the rms 

electron density fluctuation. Our result is consistent with the 

upper limit of 10* 3 4 cm~^ for AN determined by the angular broadening 

of some radio sources during occultation by the comet Arend-Roland 

(1956f) (Whitfield and Hogbom 1957). It is of interest to compare 

3 -3 

AN obtained here with the mean plasma density <N g > = 10 cm 

3 -3 

estimated from optical data by Wurm(1968). For <N g > * 10 cm and 
AN = 200 cm"\ the ratio ( AN/<N g > ) is about 0.2. 

The occultation observations provide us with valuable information 
about the plasma irregularities in the cometary tail. In order to 
determine the shape of the power spectrum of the electron density 
fluctuation, more observations with multiple channels are suggested. 
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Figure Caption 


Figure 1 . A schematic sketch of the thin screen diffraction problem. 
The random medium is confined to a thin layer of thickness D (from 
z 5 -D to z 8 0). The plane wave e 1 ^ kz_wt ) hits the "thin screen" 
from the -z direction. After passing the screen, the phase of the 
wave is randomized and is characterized by the function $(x,y). 
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